arXiv:1502.06429vl [quant-ph] 23 Feb 2015 


Quantum optical non-linearities induced by Rydberg-Rydberg interactions : 

a perturbative approach 

A. Grankin^, E. Brion^, E. Bimbard^, R. Boddeda^, I. Usmani^, A. Ourjoumtsev^, R Grangier^ 

^ Laboratoire Charles Fabry, Institut d’Optique, CNRS, 

Univ. Paris-Sud, 2 Avenue Fresnel, 91127 Palaiseau, France 
^Laboratoire Aime Cotton, CNRS, Universite Paris Sud, ENS Cachan, 91405 Orsay, France. 

In this article, we theoretically study the quantum statistical properties of the light transmitted 
through or reflected from an optical cavity, filled by an atomic medium with strong optical non¬ 
linearity induced by Rydberg-Rydberg van der Waals interactions. Atoms are driven on a two-photon 
transition from their ground state to a Rydberg level via an intermediate state by the combination 
of a weak signal field and a strong control beam. By using a perturbative approach, we get analytic 
results which remain valid in the regime of weak feeding fields, even when the intermediate state 
becomes resonant. Therefore they allow us to investigate quantitatively new features associated with 
the resonant behaviour of the system. We also propose an effective non-linear three-boson model of 
the system which, in addition to leading to the same analytic results as the original problem, sheds 
light on the physical processes at work in the system. 


I. INTRODUCTION 

Photons are considered as the best long-range quantum 
information carriers; they, however, do not directly inter¬ 
act with each other, which makes the processing of the 
information they carry problematic [l|. Standard Kerr 
dispersive non linearities obtained in non-interacting 
atomic ensembles, either in off-resonant two-level or reso¬ 
nant three-level configurations involving Electromagnet- 
ically Induced Transparency (EIT), are usually too small 
to allow for quantum non-linear optical manipulations. 
Among other techniques [ij, a possible way to enhance 
the non-linear susceptibility is to resort to a Rydberg 
level as one of the long-lived states involved in the EIT 
process dil : in such Rydberg EIT protocols, the strong 
van der Waals interactions between Rydberg atoms cre¬ 
ate a cooperative Rydberg blockade phenomenon 
where each Rydberg atom prevents the excitation of its 
neighbors inside a "blockade sphere" and deeply changes 
the EIT profile. In particular, giant dispersive non-linear 
effects were experimentally obtained in an off-resonant 
Rydberg-EIT scheme using cold rubidium atoms ^aced 
in an optical cavity [Tl| . In a previous paper @ , we 
theoretically investigated the quantum statistical prop¬ 
erties of the light generated by this scheme in the disper¬ 
sive regime, i.e. for strongly detuned intermediate state. 
We showed that, under some assumptions, the system 
effectively behaves as a large spin coupled to the cav¬ 
ity mode [l2| and we computed the steady-state second- 
order correlation function to characterize the bunched or 
antibunched emission of photons out of the cavity. 

In the present paper, we deal with the same system, 
but in a different approach. Restricting ourselves to the 
low feeding regime, we present an analytic derivation of 
the correlation function (r) for the transmitted and 
reflected light, based on the factorization of the lowest 
perturbative order of operator product averages. It is 
important to note that this derivation is valid in both 
the dispersive and resonant regimes and therefore gener¬ 


alizes our previous results. This factorization property 
is demonstrated rigorously for purely radiative damping, 
but we show also that it is approximately preserved in the 
experimentally relevant case of additional dephasing due 
to, e.g., laser frequency and intensity noise. In addition, 
we propose an effective non-linear three-boson model for 
the coupled atom-cavity system which allows us to ob¬ 
tain the same results as the (more cumbersome) exhaus¬ 
tive treatment. In the dispersive regime, this Hamilto¬ 
nian agrees with the one we obtained in the so-called 
“Rydberg-bubble approximation” it also allows us to 
investigate the dissipation at work in the resonant case. 

The paper is structured as follows. In Sec. m we re¬ 
call our setup and the assumptions we make to compute 
its dynamics. In Sec. uni we present an analytical way 
to obtain the correlation functions for the light outgoing 
from the cavity and discuss some of the numerical re¬ 
sults we obtained. In Sec. IIVI we present and discuss an 
effective three-boson model, allowing us to recover and 
generalize the previous results. Finally, we conclude in 
Sec. m by evoking open questions and perspectives of 
our work. Appendices address supplementary technical 
details which are omitted in the text for readability. 


II. THE SYSTEM 

The system we consider here is the same we dealt with 
in [H. It comprises N atoms which present a three- 
level ladder structure with a ground \g), intermediate 
|e) and Rydberg states |r) (see Figure [J). The energy 
of the atomic level \k = g,e,r) is denoted by huk (by 
convention ujg = 0) and the dipole decay rates are 7e 
(intermediate state) and 7^ (Rydberg state). The tran¬ 
sitions \g) o |e) and |e) O |r) are respectively driven 
by a weak probe field of frequency tOp and a strong 
control field of frequency Wc/. Both fields can a pri¬ 
ori be resonant or not with atomic transitions, the re¬ 
spective detunings being defined by Ag = (wp — oJe) and 
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Ar = (ujp + iOcf — uJr)- Moreover, the atoms are placed 
in an optical cavity: we shall denote by the respec¬ 

tive decay rates through the left and right mirrors (see 

Fig. [T|), with 7 c = 7 c’^^ -I- 7 !^^ The transition |g) o |e) 
is supposed in the neighborhood of a cavity resonance. 
The frequency and annihilation operator of the corre¬ 
sponding mode are denoted by ujc and a, respectively ; 
the detuning of this mode with the probe laser is defined 
by Ac = (wp — ojc) and a denotes the feeding rate of the 
cavity mode with the probe field, which is supposed real 
for simplicity. Finally, we introduce g and flc/ which 
are the single-atom coupling constant of the transition 
\g) -O' |e) with the cavity mode and the Rabi frequency 
of the control field on the transition |e) O |r), respec¬ 
tively. As represented on Figure [I] the setup allows one 
to measure the statistics of both the reflected and trans¬ 
mitted lights, i.e. g^^^ (r). 


FIG. 1. a) The setup consists of N cold atoms placed in an 
optical cavity which is fed by a weak (classical) laser beam 
of frequency ojp and a strong control laser field of frequency 
LJcf- b) The atoms present a three-level ladder structure 
{|(?), |e) , |r)}. The transitions |g) •<->• |e) and |e) -O- |r) are 
driven by the injected probe and control laser fields, respec¬ 
tively, with the respective coupling strength and Rabi fre¬ 
quency g and Qcf (see the text for the definitions of the dif¬ 
ferent detunings represented here). 


The dynamics of the full system, including the bath 
modes, are governed by the Hamiltonian derived in Ap- 
pendix|3 in the Rotating Wave Approximation. We note 
that this Hamiltonian description does not take into ac¬ 
count any additional dephasing due to, e.g., laser inten¬ 
sity or frequency fluctuations : decays and dephasing are 
therefore purely radiative. In the Markov approximation, 
the corresponding Heisenberg-Langevin equations are 


^ a = iDcO -ia-igY^ + \[ 2 ^a'b"' + 
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( 2 ) 

( 3 ) 

( 4 ) 


where and denote Langevin forces associated to the incoming fields from left and right sides and to the 

atomic operator respectively. We also introduced the complex effective detunings = (A^ -|- iy^) for k = c,e,r 
and Der = (A^ — Ae) + i ( 7 ^ + 7e)- Note that we chose to make the feeding factor a appear explicitly in Eq.([T|): in 
technical terms, it corresponds to displacing the incoming field from the coherent state ja) to the vacuum | 0 ) ; to be 
consistent with this choice, from now on, we must set (am) = 0 . 


In the next section, we show how to compute the corre- parameter a for the transmitted and reflected light, 
lation function g^'^^ (t) at the lowest order in the feeding 
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III. PERTURBATIVE CALCULATION OF g 


( 2 ) 


B. Factorization in the perturbative limit 


A. Correlation functions of the transmitted and 
reflected light. 

The second-order correlation function characterizes 

the bunched (0) > ('^))) or anti-bunched 

(0) < (5^^^ ('’’))) nature of the light transmitted or 

reflected by the cavity. For the transmitted light on 

the right side [R) of the cavity, one has by definition 

= I {R)\ iR)\ (R) {R)\ , I (R)\ {R)\^ , 

9 t v^/ — \^out ^out ^out^out j / \^out ^out j •) wnere 

is the transmitted mode field annihilation opera¬ 
tor, and all averages should be evaluated in the steady 
state. From the input-output relations (T^, one gets 

aiul + ^ and hence 

(0) = (o^o^aa) / (a^a) . 

For the reflected light on the left side (L) of the cavity, 

one gets (0) = a^l'^ , 

where is the reflected mode field annihilation oper¬ 
ator. Similarly, by using the input-output relation 


AL) 


(L) . 

aL — 1 


27. 


(i) 






for the left mirror, one gets 


= ( 27 c^^) (a^a^aa) -f 
[(a'^a'^a) - (a’^aa)] + i-^ ((a^) - (a)) -f 

7c 

d 

(4 (a^a) — (a^a^) + (aa)) + 


In the whole paper, we shall restrict ourselves to the 
low excitation regime, i.e. to low values of the feeding 
parameter a. We therefore seek ( 0 ) at the lowest non¬ 
vanishing order in a: this requires to evaluate (a^a^aa), 
and at the fourth, third and second or¬ 

ders, respectively. This task is greatly simplihed by the 
following remarkable factorization property, established 
in Appendix [B] 

(a^ (t) a = (a^ (a 

(a'l' (t2) a'*' (ti) a = (a'^' (^2) x (a 

(a’l' {12) a'*' (ti) a (ti) a (^2))*'"^^ = 

(a'l' (^2) a'l' X (o (ti) a (t2))^^^ 

where the superscript {k) denotes the order in a to which 
quantities are calculated. Therefore, for instance, for the 
transmitted light. 


5 ^^ ( 0 ) = (a) 


( 1 ) 


and we merely need to determine and Note 

that the factorization does not apply to products of the 

kind so that ^ (a) 

The mean values (a)*'^^and are readily ob¬ 

tained through taking the steady state of the first-order 
averaged Heisenberg equations Eqs. (fTTO 


= 


Dr- 


q^N 


De- 


4 B 7 


ag 


Dc {D,-^)-g^N 


aitNiul) = 27^^’ (a^a) -f ia ((a^) - (a)) -f 


27, 


(L) 


( 5 ) 

( 6 ) 
( 7 ) 


( n\ 

The second-order value (a^) is determined through 
solving the following closed system 


jA) \ — 
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(ac)^^^ 

(6c) 


g^/N , ,> (2) , a (1) 

^{aby> + -Ja)^> 

__/ac)(2) + d'/N ( 2 ) gy/N ( 2 ) _a_ 

2(7^, + 7^e)^ ^ {D, + D,)^ ‘ ^ {D,^De)^ ' {D, + De) 


g'/N (2) a 

[D, + i?,) ^ ^ (Z?, + Dr) 


(c)(^) + 


(66)^ 

: -I- 


(&) 


( 1 ) 


2 (-Dc + Dr) 


(ab)^ 


n 


cf 


2 (-Dg + ^r) 


(cc)*-^^ + 


g'/N 

{De + Dr) 


(ac) 


( 2 ) 




c/ 


2 (z?e + -Dr) 


( 66 ) 


( 2 ) 


^cfgy/^ (2) ^^fg^ ^ (2) 

^ (ac)(2) + -^%—K (a6)(") 


( 8 ) 

( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 

( 13 ) 


deduced from Eqs. (ITEl) under the assumption of an homogeneous atomic medium, whose consequences are detailed 
in Appendix [C] In this system, we introduced the collective atomic operators 
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We note that the first-order mean values (a)^^\ {b)^^'^ 
and (c)*'^^ which appear in Eqs. (Eiiaiini), respectively, 
have been computed in Eqs. ©[mD. The K coefficient 
is approximately given by (see Appendix [C] for details) 


K R 

where 


' jg + £) _ Zs/ 


r)2 

n _ ifb 
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1 - 


ypiir^ 
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-C. 


Dr - Dlfl ( 4 (De + Dr) - 


( 14 ) 


( 15 ) 


will be interpreted as the Rydberg bubble volume in 
the dispersive regime in the next section. Though it 
is too cumbersome to be reproduced here, the solution 

^ 2 ) 

for {a^') is simply obtained by matrix inversion, and 

the calculation of (0) and (0) can be straightfor¬ 
wardly programmed, e.g. in Mathematica. 

As it has been the case for g[^r (0)) Ih® calculation 

(‘2'\ 

of the time-dependent correlation function ^ (r) = 

(al (t) al (t -I- r) a (t -I- r) a (t)) / (ala)^ is greatly simpli¬ 
fied by the factorization property derived in Appendix 
m since we simply need to determine the quantity 
(a (t -I- r) o (t)). From Eqs. (II]l 4 |), one easily deduces the 
following differential system, at the lowest order in a. 


, / {a{t + T)a{t)) 

f (6(t + r)a(t)) 
^ \ (c {t + T)a (t)) 


—ia (a) 



-b 


/ -Dr gVN 0 \ / {a (t + T) a {t)) \ 

-i g^/N -Dr ^ W {b{t + T)a{t)) 

\ 0 -Dr ) \ h + ) 


which, together with the initial condition 


/ {a{t + T)a{t)) 
{b{t + T)a{t)) 

V {c{t + T)a{t)) 



( \ 

- 

(ba)^^^ 

f t=Q 

\ (ca)^^^ / 


calculated above, allows us to determine (a (t -I- r) a (t)). 
Again, though involved, the expressions are straightfor¬ 
ward to obtain and program. 


C. Application to an experimental case. 

1. Dispersive regime. 


Let us now provide some numerical results obtained in 
the perturbative approach described above. We first in¬ 
vestigate the dispersive non-resonant regime, addressed 
in our previous work To be specific, we consider the 
same system, namely an ensemble of ®^Rb atoms, whose 


state space is restricted to the levels \g) = 


5si;F = 2), 


e) = 5p|;E = 3^ and |r) = 


95d=;E = 4 


The re¬ 


spective radiative decay rates are 7 e = 27r x 3 MHz and 
7 r = 27 r X 0.03 MHz, the cavity decay rate is 7 c = 27 r x 1 
MHz, the volume of the sample is H = 407r x 15^ Itni^, 
the sample density Uat = 0.4 /rm“^, and the cooperativ- 
ity C = g‘^N/{2jelc) = 1000. 

The other parameters take the same values as in in 
units of 7 e, the control laser Rabi frequency is Qcf = 10, 
the detuning of the intermediate level is Ac = —35, the 
detuning of the Rydberg level is A,. = 0.4, the cavity 
feeding rate is a = 0.01, and the Van der Waals coeffi¬ 
cient is Cq = —8.83 X 10®7e p,m®. For these parameters, 
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FIG. 2. Second-order correlation function (0) for the 
transmitted light in the dispersive regime considered in @ 
as a function of the renormalized cavity detuning Ocj^e = 
/ 7 e where Ac°^ is the detuning of the linear cav¬ 
ity. The shape of the plot is in good qualitative agreement 
with the results of the previous model. Inset : the same plot 
in logarithmic scale (bunching and antibunching peaks are 
more clearly visible). 


the maximal average number of photons in the cavity 
is obtained for the cavity detuning = —6.15206 7 e 
which is taken as a reference. 

Let us note however that in real experimental con¬ 
ditions, the atoms undergo not only radiative damp¬ 
ing, but are also subject to extra dephasing on the 
Rydberg-ground state transition, due to laser frequency 
and intensity noise. This additional dephasing cannot 
be modeled in the Hamiltonian formalism presented in 
Appendix El and thus the demonstration given in Ap¬ 
pendix m for the factorization of mean values does not 
apply any more. However, since the radiative coher¬ 
ence damping is jr ~ 0 . 017 e, the additional damping is 
'yd ~ 0.15 7 e, and the total number of atoms in the sam¬ 
ple is N Ki 10"^, the experimental parameters satisfy the 
condition 7 r 7 d <C N"/r. Under these circumstances, 
it is shown in Appendix |D] that the factorization remains 
valid, provided that the coherence radiative damping 'yr 
is replaced by the dephasing rate in the equations. 

Under these conditions. Figure [H shows the second- 
order correlation function (0) as a function of the 

reduced cavity detuning 6c = (Ac — A^^^j/yg, to be com¬ 
pared with Fig. 2 a) in Q. The two plots are in good 
qualitative agreement, but the position of the bunching 
peak is shifted from 6c ~ —5 to 6c ~ —3.5, for the same 
parameters. This basically originates from the definition 
of 14 in @ , differing from the present one by a factor v^. 


before. As a new set of paramenters, we take Ac = Ac = 
Ac = 0, and we assume that yc^^ <C 7c^^. We also 
choose a higher principal number n = 100 for the Ry¬ 
dberg level, for which y^ = O.lyg. In addition, we fix 
y^ = 0.3ye, C = « 30 and U = 507r x 20 x 20/im3. 

In this regime, 14 ~ enhanced, therefore 

magnified non-linear effects are expected. 


i^out^ut^out^out) > 



FIG. 3. Resonant case A^ = Ae = A^ = 0. The quanti¬ 
ties (thick line) and (thin 

line), renormalized by the intensity of the incoming light, are 
represented as functions of the normalized control field Rabi 
frequency Qcf /'Je- For flcf = 2-\/yeyr (2(7 — 1) « 5ye, photon 
pairs are reflected, i.e. / 0, while single 

photons are absorbed, i.e. ~ 0- 


As can be seen on Figure |31 there exists a value for 
which single photons are mostly absorbed = 

0, while pairs are reflected ^ 0: this 

value can be computed and is found to be 


He/ = (2C - 1) = 2ye'\/6 « bye 


2. Resonant case 

After checking that the present work confirms our pre¬ 
vious results, obtained in the dispersive regime, let us 
consider the resonant case, which could not be treated 


On the contrary, in a slightly detuned case, i.e. for Ac = 
— 2ye and A^ = — O.lye, the other parameters remaining 
the same, one observes that around Hg/ Ri llyg pairs are 

absorbed = 0 while single photons 

are reflected 0 (see Fig. i]). 
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FIG. 4. Slightly detuned case Ac = 0, Ae = —27e, = 

-0.l7e. The quantities (thick line) and 

(thin line), renormalized by the intensity of the 
incoming light, are represented as functions of the normalized 
control field Rabi frequency flc//7e- For Q^f ~ ll7e, photon 
pairs are absorbed, i.e. = 0, while single 

photons are reflected, i.e. 7^ 0- 


These new features are specific of the near-resonant 
regime, and were not present in our previous work. They 
may be interpreted as different impedance matching con¬ 
ditions for single photons and for pairs, leading to very 
large non-linear losses, acting at the single photon level. 

To conclude this section, we described how to obtain 
the exact and analytic expression of the correlation func¬ 
tion in the low excitation regime, valid not only in the 
dispersive regime but even in the resonant case. Though 
exact and computable, the expressions we get are too 
cumbersome to be displayed here and do not easily lend 
themselves to physical interpretation. In the next sec¬ 
tion, we introduce an effective non-linear three-boson 
model which allows us to derive the same results to the 
lowest order, and has also the advantage of being physi¬ 
cally more transparent. 


tonian is given by 


H = — Actt^a -I- a (a -I- a^) — Agb^b 
— Arc'c + g'/N {ah"^ -|- &^a) 

+ -I- ^c^c^cc 

We moreover assume that the c-boson is coupled to a 
non-linear bath whose action on the system is repre¬ 
sented by the following non-linear dissipation operator, 
acting on the density matrix p of the system 


'Dip] = -^ {2ccpc^c^ — c^c^ccp — pc^c^cc} 


Here, all parameters, in particular and are assumed 
real. From the full Liouville-von Neumann equation of 
the system dtp = — ^ [H,p] D[p] one readily derives 
the following Bloch equations 


^ (a) = iDc (a) -ia- igy/N {b) 

^ (6) = iDe {b) - ig'/N (a) - i^ (c) 
^ (c) = iDr (c) - i^ {b) - in (c+cc) 


IV. EFFECTIVE NON-LINEAR THREE-BOSON 
MODEL 

A. Non-linear absorption and dispersion in the 
quantum regime. 

We consider a system of three bosons of respective an¬ 
nihilation operators a, b and c, whose non-linear Hamil- 


where we introduced the notation k = Kr — iKi. From 
this set of equations, one gets the same steady state value 
as in Eq. At the second order in a, the set 

of equations for two-operator steady-state averages is de¬ 
rived in the same way (here we omit superscripts for 
simplicity) 
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(ac) 
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{be) 

(cc) 


-ft- (“'■> + IT <“> 




cf 


2 {Dc + De) 

gy/N 


/ \ , 5^^ / \ , , 


{Dc + Dr) 


{he) 


{Dc + De) 
a 


^cf ,j,„\ , gVN 




( 7 ?c + Dr) 
{ab) 


(c) 


{Dc + De) 

D,cf 


{Dc + De) 


(b) 


2 {Dc + Dr) 


{ab) 


a 


cf 


2 (-Dfi + Dr) 

D,cf 


De 

/ \ , dVN 

(cc) + —-—— (ac) + 




c/ 


{De + Dr) 


2 (Z^r 2 ) 


(6c) 


2 {De + Dr) 


( 66 ) 


which agrees with Eqs. (l8l fT^ but for the last equation. 
If, however, we eliminate(6c) and (66) from the last three 
equations, one obtains 


(cc) = 




n n_ sL ] _£l 

-t- 4£) I 4 


^ (ai.) 


which can be identified with Eq. (ITSl) provided that 

1 


K = 


(Dr-^) (^Dr+De-^)-^ 


which, upon recalling Eq. (HI, yields 


K = 2 


Vb 


V-Vb 




4 ( Dr + De — 


4D^J 


- Dr 


the Rydberg bubble approximation @ and therefore con¬ 
firms its validity: we observe a shift due to the non-linear 
dispersive behavior of the c-boson, but no non-linear ab¬ 
sorption since the intermediate level is too far detuned. 
Moreover, in the bubble picture, the parameter Nb was 
interpreted as the number of Rydberg bubbles the sam¬ 
ple may accommodate; as suggested above, this allows to 
interpret 14 as the bubble volume. 

If we now go to the opposite regime, i.e. the resonant 
case for which Ag = = 0 , 7 e 7 r and 7 ^, 

we obtain I 4 ~ ^ (I — i) -\J and therefore the non¬ 
linearity parameters are 

27r2 _ 

Hr = —Hi ~ \rle\C^ 

We now have both dispersion and absorption. From the 
expression of it is clear that absorption results from 
an interplay of the spontaneous emission from the inter¬ 
mediate state and the Rydberg-Rydberg interactions. 


V. CONCLUSION 


We obtain thus the analytic expressions of the parame¬ 
ters Kr and Kj, respectively characterizing the non-linear 
dispersion and absorption of the c-boson, which make our 
model system precisely reproduce the results of the orig¬ 
inal problem in the steady state and in the lowest order 
of the feeding parameter a. 

B. Discussion. 

Let us now investigate the physical content of the pre¬ 
vious model by considering two limiting cases. 

In the dispersive regime addressed in our previous 

work \De,r\ ^ f^c/, whence 14 ~ ~ 

^-iid Ki ~ 0, where we introduced Nb = 

This result agrees with what we previously obtained in 


In this article, we have studied the strong quantum 
optical non-linearities induced by Rydberg-Rydberg van 
der Waals interactions in an atomic medium. We pro¬ 
vided a new perturbative treatment of the problem, based 
on the factorization of the lowest perturbative order of 
operator product averages. Though being purely ra¬ 
diative damping, this factorization property is approx¬ 
imately preserved in the presence of to, e.g., laser fre¬ 
quency and intensity noise, as it is the case in our exper¬ 
imental setup. Our perturbative calculations enabled us 
to recover and extend our previous results: we could vali¬ 
date the approach based on the Rydberg bubble picture, 
as well as investigate the resonant, absorptive, regime. 
In particular, our numerical simulations showed that 
strong Rydberg-induced non-linearities led to different 
impedance matching conditions for single photons and 
photon pairs. 
























Moreover we proposed an effective model which leads 
to the same results as the full calculation at the lowest or¬ 
der in the feeding parameter; this model also sheds some 
light on the origin of the dispersion and absorption, as 
well as makes a bridge between the Rydberg bubble and 
perturbative approaches. In the future, we shall first try 
and take advantage of our understanding of the system 
to investigate regimes of parameters for which a photonic 
gate can be implemented. On the other hand, we also 
plan to apply other methods, inspired from many-body 
physics to the problem, in order to recover and further 
extend our results. 
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Appendix A: The full Hamiltonian in the Rotating 
Wave Approximation 

The full Hamiltonian of the system can be written un¬ 
der the form 

H = Hat + Hcav + Hbath + Vat — cav ^cav— bath — bath 

N N 

Hat = fy^e ^ + Klcf COS {uJcft) ^ + 

n=l n—1 

N N 

n—1 m<.n—l 

Hcav = huJctt^a 

Hbath = J doj huj (b^b + ch + d^d') 

N 

Vat-cav = '^hg{a + a^) 
n—1 

Vcav—bath — j duJ figb {p 4“ (a U^) 

Vat-bath = J ^9c (c -k c'l') -k -k 

J duj hgd {d + d'^) 

where Uap = |a) (/3|, Huj a is the energy of the atomic 
level I a) for a = e,r (with the convention Wg =0), and 
Kmn = Oe/ ||#n — r„||® denotes the van der Waals inter¬ 
action between atoms in the Rydberg level - when atoms 
are in the ground or intermediate states, their interac¬ 
tions are neglected. The operators b{uj), c (w) and d{uj) 
denoted simply as b, c, d, are bath operators coupled 
to the cavity and atomic operators with the respective 
coupling strengths gb {to), gc (to) and gd{to). 

We switch to the rotating frame defined by |#) —>■ 


= exp {—^Hq) where 

N 

Ho = hjOpO^a -k ^ {^pCT^JP -k h (tOp + UOcf) + 

n—1 

and perform the Rotating Wave Approximation to get 
the new Hamiltonian 

H = Hat + Hcav + Hbath + Vat-cav + Vat — bath ^cav—bath 

N N 

Hat = -hAcY,<^i:'^ + E + 


n—1 m<.n—l 


N fcO ^ 

-KA.j: 4 k+#^i:( 4 ;>+dk) 


n—1 

Hcav — ^ 


n—1 


Hbath ~ / doo hto {b{to + Wp))^ b{u} + uop) + 

f ^ 

dto fko E (Ctx (w “1“ tOp ~h UJcf')') ^n “t“ ^p ~\~ ^cf') “t“ 

n^l 

AT 

diu huJ ^ {dn (oj + dn {uj + 0;^) 

n^l 

N 

Va-c ~ E 

n—1 

^cav — bath ~ / dto hgb (to) 1^6 (to) -k (b (w))^ aj 

w r r 

Vat-bath j duo hgc (uo) Cn (uo) -k (c„ (w))''' 

^ r r 

/ dto Hgd (uo) dn (to) aif + (dn (uj))^ 


n—1 


with the detunings Ac = (tOp — uoc), Ag = (tOp — toc), and 
Ac = (tOp + tOcf — tOr). It is important to note that the 
evolution under the Hamiltonian H conserves the number 
of excitations. 


Appendix B: Factorization of correlation functions. 

We suppose that the bath interacting with the cavity is 
initially in the following continuous-mode coherent state 
(incoming quasi-classical field) 

|a) |0) 

where f la (t)P dt = (n) and 6 ^ = / f dtoa (to) b^ (uo) 

vO> 

is a superposition of bath mode creation operators b^ (to) 
M- Note that with this definition, is a bosonic opera¬ 
tor, i.e. [ 6 a, 6 y = 1. The atoms and cavity field are ini¬ 
tially in their ground state denoted by |G) = |g ... 5 )( 8 )| 0 ). 
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Let us consider, for instance, the quantity 
(a,G\a^ {ti)a^ {t 2 )a{t 2 )a{ti)\G,a), for t 2 > ti, 
where jG, a) denotes the initial state of the whole 
system {atoms+cavity+baths}, the baths coupled to 
the atoms are supposed empty and their state is not 
explicitly written, 

(a, G|a'^ (ti) a'^ (^2) a (^2) a (G) |G, a) (Bl) 

k + l 

Expanding this expression with respect to |q;| (which is 
equivalent to expanding in the number of excitations 
present in the system), one finds that the lowest non¬ 
vanishing contribution is the fourth order term k = I = 2. 
For the system considered the identity operator can be 
represented in the following way I = where li = 

i 

k*) (ft I til® identity operators on each degree of 
freedom of the system, and \qi)^s denote q -th basis vec¬ 
tor of i -th degree of freedom. Inserting this identity 
operator between a'*' (t2) and a (t2) of the quantity (IBII) 
yields: 

k + l 

e-(-)^!^{k,G\aHh)aHt2) (B2) 

<S> X! ki) (ftl I a (^2) a (G) |G, 1) 

i q ) 

For the lowest non-vanishing term A: = 2, / = 2: 

a {t 2 ) a ih) |G, 2) = |G, 2 (G)) 

where |G, 2(ti)) = e“‘-R^|G, 2) (note that this state 
can contain excited atoms and/or cavity photons). 
The state a|G, 2(ti)) can at most contain one ex¬ 
citation, and so can the state e’’^ (ti)) 

due to the conservation of excitation number. Hence 
e‘^ 2 (H)) can only have component on 

|G, 0). Finally the fourth order expression of (|B2I) reads: 

e-M (2, G|at (G) (^2) a (^2)« (G) |G, 2) 

= e-w(!|_|(2,G|at (G) (t2)|G,0)|" 

= (a, G|al (H) a) { 12 ) |G, O)^ (G, 0|a (^2) a (ti) |G, a)^ 


where we used that e ^ <^2, G|al (ti) (^2) |G, O) 

and (0, G|a (^2) a (ti) |G, 2) are equal to 

the second order expansion in |a| of quantities 
(a, G|al (<i) (^2) |G, a) and (a, G|al (ti) (t2) |G, a) 

respectively, which we denoted by {■■■) 2 - 

Thus to compute (a, G|al (G) (G) a (G) a (G) |G, a) 

in the lowest order it is enough to calculate 
(a (G) a (ti)) = (a, G|a (G)« (G) |G, a). 

The same argument holds for more general mean 
values such as 

{a, G\ al (G) o' (G) (tp) a (tp+i )... 

... a (G-i-g-i) a {tp+q) |G, 

and in particular 

(al (t) a = (al (t))^^^ (a 

(al (G) a} (G) a (ti))*'^^ = {a^ (G) (G))*'^^ 

X (a (G))^^^ 

(al (G) (G) a (G) a (G))^"^^ = (al (G) o' (ti))^^^ 

X (a (G)a(G))*'^^ 


Appendix C: Calculation of 


The system of equations for the same-time 2-operator 
products in the second order in a is readily derived from 
Heisenberg-Langevin equations. For notational conve¬ 
nience here, we do not explicitly write superscripts 
nor the time since we only dealt with same-time mean val¬ 
ues : hence (aa) should be understood as (a (t) a (<))^^^ 

and (^cTge'^ as ^fJge (t)^ ■ We thus find 
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_d 

dt 



d 

dt 



= 2Dc {aa) — 2ig'^^ — 2ia (a) 

i 



3 



Assuming that the medium is homogeneous, i.e. that 
for all and (aa^g}^ = 

I 


r 


^acTgi^^, in the steady state this system yields 


(»d 2 ) = 2 


Se 




a 


_ 

{D, + Dr)^ • 


3 

iD, + Dr) 


gr 


\^ge ^ge j 
ge ^gr j 
\^gr ^gr / 


■'c ' ' 

i 

nyW *!“"”) + JD^) <““> + (zyfiy ^ 
(D, + D,) S (di’d? 

-- /cr(^'V(*A - 

2 (lie + Dr) \ 

Qi. / /.^ ,.s\ 


a 

(Dc + 


De) 






5 


Note that the first-order values (a) = (a)^^\ — 




( 1 ) 


Ji) 


( 1 ) 


have been determined 


through solving the first-order steady state system, see 
Eqs. ([5]l7| in the main text. 


we get 


a(^V«) = 

^ gr ^ gr ' 


ntg 


Dr - {Df, + Dr) — 

^l9 


ii. 

4De 


■} 


aa, 


{i) 


^D,[[Dr-^) [[D, + Dr)-§-]-^] 


0 ( 7 , 


(0 


Summing the above equations over atom numbers (i, j) 
yields a system on averages of the collective operators 
b = c = ^J2i ^gr and field operator 

o, which is almost closed but for the last equation which 
will now be considered and approximated. Eliminating 

cr^e ag^)\ and (crgiVge / from the last three equations 


We now sum over i and j indices and divide by N this 
equation to get 
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where we introduced the coefhcient 




D, 


Dr - 


A) 


4 


Making the approximation that Ki does not depend on 
i, i.e. Ki m K, we get: 


K =—Y\ 

N ^ 


47r 

47r 


[De + Dr-A) (Dr-^)-^ 


i?3 


/ 


A. + Dr-A) {Dr-^)- 


i?3 


/' 


( 


De+Dr- 


JZL'i (D — 

AD^J 2rB) 4 


£i 

4 

dr 


-dr 


For large values of R, K does not depend on the geometry 


(cc) 


2 


K (ac) + 


^IgVN 


K (ab) 


To estimate K we consider that the sample is a sphere of 
radius R 

I 


(D + D _ D — ^ 

A- J-Jr 4JJ^ J Ayr 4 


/ 


\ 

1 - 


Ce 

3H 

A 

p 

V 


\De+Dr-^j y 


Finally the desired closed system is 


(aa) 

(ab) 


(ac) 

m 

{be) 


(cc) 


9'^ / r\ , “ / \ 

-57 <“*’> + d: <“> 


Q.h 


2 {Dc + Dg) 

gy/N 


/ \ , dy/N gy/N 

(“c) + /n , n ^ (aa) + /n , n N 


{Dc + Df.) 


{Dc + Dr) 


{be) + 


{Dc + Dr) 


{c) + 


{Dc + Dc) 

rif, 


(He + Dc) 


(b) 


2 (He + Dr) 


{ab) 


> gy/N 

Jd. + -ST 


(ac) + ^ (ac) + 


rif, 


2 {Dc + Dr) 
nbgVN 


{Dc + He) 


2 (He + Dr) 


(56) 


iF (ac) + 


nlgy/N 


iDc 


K {ab) 


(Cl) 


which allows to determine {aa). The analytical solu¬ 
tion is too cumbersome to be displayed in this paper but 
can be readily obtained by matrix inversion. 

Appendix D: Factorization in the presence of extra 
dephasing 

In this appendix, we show in which conditions the 
factorization of field operator products described in Ap¬ 
pendix |B] remains valid in the presence of extra dephas¬ 
ing due to laser frequency and intensity noise. Such de¬ 
phasing is correctly accounted for by adding the term 
—+ Hg? in the Heisenberg-Langevin equation Eq. 
([5]) on (Tgr\ where Hgf^ is an extra Langevin force and 


7d ~ 0.15 X 7e, 7r- ~ 0.01 X 7e and 7e = 27r x 3 MHz in 
the experimental setup. 

In the absence of interatomic interactions, because 
laser and cavity fields address the atoms symmetrically, 
the ensemble evolves in the subspace of symmetric states. 
The atomic system essentially remains in this subspace, 
even when the interactions are taken into account, if the 
number of Rydberg excitations in the sample is much less 
than the total number of Rydberg bubbles the ensemble 
can accomodate for. Such symmetric superpositions ac¬ 
tually not only contain “allowed” components {i.e. with 
Rydberg atoms further than a Rydberg bubble radius 
apart) but also “forbidden” components (with Rydberg 
atoms closer than a Rydberg bubble radius). Their num¬ 
ber is, however, very small compared to that of “allowed” 
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confgurations and they will therefore only slightly alter 
the outcome of dissipative dynamics of the system. 

Under these assumptions, let us show in which con¬ 
ditions the mean value (c^c) factorizes at lowest order. 
Focusing on the dissipative part of Bloch equations for 
(Tgr and CTrr (note that for the latter, there is no extra 
dephasing) we get 


dt 


dt 



d 



yU) 


d,i¥=j 



- {ir + Id) ) 

-2(7r +7<i) {^rg^^ir) 

-2-fr <^crW ^ 


and recalling that c = we get (c'^c) = 

E* ) + ;^ Ei/J {^rg ^ whencc, for a short 
time interval 


d{c^c) 


dt 


dt 


N dd ' N dd ^ ^gr / 




= ^ (7r + Id) 


N 


N 




= -^ ) +1 (>+)+ 


- ^ (7r + 7d) '^{(^rg^gr) 


N 

27d 




d- ^ E^^") “^(7r- + 7d)(c^c) 




When there are Ur Rydberg excitations in the sample, 
with Tir «C Nb <C N (Nh is the maximum number of 
Rydberg excitations the sample can contain), one has 

(c^c) « Ei (cTrr \ « nr whencc 



-2 


7r +7d 



(c^c) 


and for 7r <C 7d "C A^7r 

^ -27d(ctc) 


so, from the point of view of c^c, everything works as if 
the system was radiatively damped with the rate 7^. In 
the same conditions, we moreover have 



-Id (c) 


and again, from the point of view of c, everything works 
as if the system was radiatively damped with the rate 7^. 
Moreover, since all other dynamical equations (for popu¬ 
lation, coherence and field operator mean values) remain 
formally the same as in the purely radiative damping, the 
factorization procedure remains valid for (a^a) provided 
that 7r 7d 'C fV7r and the radiative coherence decay 
7r is effectively replaced by the dephasing decay rate ■jd- 

This result can also be extended to higher order quan¬ 
tities ((a'l')’"aP). 


[1] D. E. Chang, V. Vuletic and M. D. Lukin, Nature Pho¬ 
tonics 8, 685 (2014). 

[2] J. D. Pritchard, D. Maxwell, A. Gauguet, K. J. Weath- 
erill, M. P. A. Jones, and C. S. Adams, Phys. Rev. Lett. 
105 , 193603 (2010). 

[3] Y. O. Dudin and A. Kuzmich, Science 336, 887 (2012). 

[4] T. Peyronel, O. Firstenberg, Q. Y. Liang, S. Hofferberth, 
A. V. Gorshkov, T. Pohl, M. D. Lukin and V. Vuletic, 
Nature 488, 57 (2012). 

[5] D. Maxwell, D. J. Szwer, D. P. Barato, H. Busche, J. D. 
Pritchard, A. Gauguet, K. J. Weatherill, M. P. A. Jones, 
and G. S. Adams, Phys. Rev. Lett. 110 , 103001 (2013). 

[6] A. Grankin, E. Brion, E. Bimbard, R. Boddeda, 1. Us- 
mani, A. Ourjoumtsev and P. Grangier, NJP 16, 043020 
(2014). 

[7] M. D. Lukin, M. Fleischhauer, R. Cote, L. M. Duan, D. 
Jaksch, J. 1. Cirac, and P. Zoller, Phys. Rev. Lett. 87, 


037901 (2001). 

[8] M. Saffman, T. G. Walker, and K. Moelmer, Rev. Mod. 
Phys. 82, 2313 (2010). 

[9] D. Gomparat and P. Pillet, J. Opt. Soc. Am. B 27, A208 

( 2010 ). 

[10] V. Parigi, E. Bimbard, J. Stanojevic, A. J. Hilliard, 
F. Nogrette, R. Tualle-Brouri, A. Ourjoumtsev, and P. 
Grangier, Phys. Rev. Lett. 109 , 233602 (2012). 

[11] J. Stanojevic, V. Parigi, E. Bimbard, A. Ourjoumtsev, 
and P. Grangier, Phys. Rev. A 88, 053845 (2013). 

[12] C. Guerlin, E. Brion, T. Esslinger, and K. Moelmer, Cav¬ 
ity quantum electrodynamics with a Rydberg-blocked 
atomic ensemble Phys. Rev. A 82, 053832 (2010). 

[13] D. F. Walls and G. J. Milburn, Quantum Optics 2nd ed., 
Springer-Verlag Berlin Heidelberg (2008). 

[14] R. Loudon, The Quantum Theory of Light, Oxford Uni¬ 
versity Press, 3rd edition (2000). 



